Adiabatic response for Lindblad dynamics 



J.E. Avron 

Department of Physics, Technion, 32000 Haifa, Israel 

M. Fraas, G.M. Graf 

Theoretische Physik, ETH Zurich, 8093 Zurich, Switzerland 

February 28, 2012 
Abstract 

We study the adiabatic response of open systems governed by Lind- 
blad evolutions. In such systems, there is an ambiguity in the assign- 
ment of observables to fluxes (rates) such as velocities and currents. 
For the appropriate notion of flux, the formulas for the transport coef- 
ficients arc simple and explicit and arc governed by the parallel trans- 
port on the manifold of instantaneous stationary states. Among our 
results we show that the response coefficients of open systems, whose 
stationary states are projections, is given by the adiabatic curvature. 

1 Introduction 

We are interested in extending the theory of adiabatic response of quantum 
systems undergoing unitary evolution [10, 23] to open (quantum) systems 
governed by Lindblad evolutions. In particular, we are interested in a geo- 
metric interpretation of the response coefficients. 

In open systems there is usually some choice in setting the boundary 
between the system and the bath. Setting the boundary fixes the tensor 
product structure Hs 8) Hb- Choosing a boundary still leaves a residual 
ambiguity in observables. For example, given a joint Hamiltonian H of the 
system and the bath, there is no unique way of assigning to H an observable 
of the form Hg 1 describing the energy of the system alone. (For example, 
the interaction of an atom with the photonic vacuum has two effects on the 



atom: It leads to decay and to "Lamb shift" of the energy levels. One can 
choose whether to incorporate the Lamb shift in the energy of the atom or in 
its interaction with the bath.) A further aspect of the ambiguity arises when 
considering the flux (rate) of an observable X. Any assignment X Xs of 
system observables Xs to joint observables X is incompatible with dynamics, 
if the bath and the system interact. In fact it is generally impossible to satisfy 
both requirements ® 1 h-> and X = i[H, X] t— )■ = i[Hs,Xs]. 

The ambiguity in fluxes is physical and plays a key role in this work. 
Consider, for example, damped harmonic motion. By Newton, the flux of 
the momentum is the total force. This force is related to two other forces in 
this problem: 

Momentum flux = Spring force + Friction force. 

The momentum flux can be determined from the trajectory of the particle; 
The spring force from the force acting on the spring anchor and the friction 
from the momentum transfer to the bath. All these forces have physical 
significance and are associated with different measurements. In this work we 
shall focus on observables that are the analog of the momentum flux. 

A time honored strategy to study open systems is to start from a complete 
Hamiltonian description of the system and the bath [21, 15, 22, 2]. This 
approach comes at the price of analytical and computational difliculties. Here 
we choose the complementary "effective" Markovian approach [3, 14] which 
is often analytically and computationally simpler but is not universally valid. 
More precisely, we study open systems in the Lindblad framework [14, 13, 17]. 

The Lindblad operator, denoted C, is made of a self-adjoint H represent- 
ing the "energy" of the system and a collection of operators, {Lq}, represent- 
ing the coupling to the bath. The notion of "energy" is, as we have noted, 
ambiguous and this is manifested in the non-uniqueness of {if, Fq,}. 

We shall call a choice {H, Fq,} a gauge. Different gauges generate the 
same dynamics. 

We shall consider parametrized Lindbladians, C^, where the (classical) 
parameters (p E Ai are viewed as controls. This means that {H{(j)),Ta{4>)} 
are functions of the controls. M. is the control space, see Fig. 1. 

The main focus of this work is the development of an adiabatic^ response 
theory for the fluxes; Namely, observables of the form X = £^(X). This no- 

^The notion of adiabaticity is contingent on a gap condition, Assumption 2 below. 
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tion is gauge invariant (independent of the choice {H, Fq,}). It turns out that 
the adiabatic response of such observables has several simphfying features. 

A different perspective on the choice of observables comes from a gauge 
invariant formulation of the principle of virtual work. For isolated systems 
the principle of virtual work assigns the observable d^H with the variation 
of the /i-th control (j)^. Since 6H is gauge dependent in the Lindblad setting, 
formulating a gauge invariant notion of the principle of virtual work requires 
the joint variation {6H, 5Fq,}. 




Figure 1: The surface represents the control space M. and the closed curve 
represents a closed path in the space of controls. The cross denotes a puncture 
in Ai where the adiabatic theory fails and the manifold of stationary states 
is singular. 

Consider a path in control space which is traversed adiabatically, Fig. 1. 
It is a feature of adiabatic evolutions [4] that stationary states evolve by 
parallel transport within the manifold of (instantaneous) stationary states. 
Often this evolution fully determines the response and transport coefficients 
for fluxes. The transport formulas can be viewed as an analog of Kubo 
formulas for linear response, emphasizing the geometric flavor. 

Although parallel transport captures the geometric aspects of the evolu- 
tion, as a practical method of calculation of transport coefficients, it suffers 
from its reliance on solving differential equations. There are, however, two 
important families of Lindbladians where parallel transport can be deter- 
mined algebraically, without solving any differential equations. This is the 
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case for (generic) Lindbladians where the (instantaneous) stationary state 
is unique. It is also the case for a second family of Lindbladians, namely, 
dephasing Lindbladians. These share the stationary states with the Hamil- 
tonian evolution [4]. In both cases the transport formulas, the analogs of 
Kubo formulas, are both geometric and explicit. 

In general, the response coefficients of open and closed systems are differ- 
ent. One would like to identify those transport coefficients that are immune 
to certain mechanism of decoherence and dephasing. For observables of the 
form £*{X) the response coefficients depend on the manifold of stationary 
states (but not on the underlying dynamics). Immunity then follows when- 
ever the stationary states are unaffected by decoherence and dephasing. This 
is the case for two physically interesting families of Lindbladians: Dephasing 
Lindbladians and Lindbladian which allow for decay to the (Hamiltonian) 
ground state. 

2 Lindbladians 

The Lindblad (super)^ operator [14, 13] is given formally by^ 

c{p)^-i[H,p] + v{p), p(p) = 5^2r«pr;-r;r,p-pr;r„, (i) 

a 

where the state p is trace class. The Hamiltonian part H is self-adjoint (and 
local). Ta are essentially arbitrary. Models describing exchange of energy 
involve non-self- adjoint while models of measurement involve Fq which 
are spectral projections (non-local in general). 

jC is the generator of state and trace preserving contractions. When T) = 
the evolution is unitary. To avoid technical difficulties with unbounded 
operators we shall assume: 

Assumption 1. 

• Weak version: H and F„ hounded 

• Strong version: H and Fq, finite dimensional. 

^ Super operators will be denoted by script characters. 
•^The normalization differs by factor 2 from that of [19]. 
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However, we shall also occasionally consider standard physical examples 
with unbounded operators. We shall study the time evolution of the state p: 

P-m, (p^l). (2) 

2.1 Gauge transformations 

C does not determine {H, Fq}. In fact, C is invariant under the joint variation 
[13] 

5H^el-iY,{9l^a-9aK). ^^a^9a. (^aGC, eGK). (3) 

a 

Moreover, V and lAT represent the same Lindbladian when lA is unitary in 
the sense that 

(wr)„ = 5^w,^r^, u-^^u*. (4) 

Wc shall refer to the freedom in {H, r^} as gauge freedom. The observable 
iJ, which one would like to interpret as the energy of the system, is therefore 
ambiguous a priori. However, in the examples we shall consider, one choice 
is singled out naturally. 

2.2 Lindbladians with a unique stationary state 

It is convenient to introduce a notation that distinguishes stationary states 
a from general states p. An important class of Lindbladians are those that 
have a unique stationary state cr with C{(7) — (with Trcr = 1). This is the 
generic situation in the finite dimensional case. The (super) projections V 
on the stationary state is given by 

V{p)^al^p, (IV(j = l). (5) 

Evidently V is trace preserving and V"^ — V (since Trcr = 1). It is not 
orthogonal, not even formally. In fact, the dual projection V* , which acts 
naturally on observables, is given by a different expression, 

V*{X)^1-T^{X(t). (6) 

A basic identity we shall need is 

£P = = VC. (7) 
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The first equality is evident. The second follows from 

P£p = aTr(£(p)) =0. 
We shall denote by Q the complementary projection 

V + Q^l. (8) 

Evidently, 

CQ^C^ QC. (9) 
This will play a role in the sequel. 

2.3 Dephasing Lindbladians 

Dephasing Lindbladians are intermediate between the unitary family and 
Lindbladians with a unique stationary state. They are characterized by = 
TaiH) for some functions Fq. In particular, 

CiP)=0 = i[H,P], (10) 

where P is a spectral projection for H. 

In the finite dimensional case all spectral projections Pj are finite di- 
mensional and dephasing Lindbladians share the stationary states with the 
Hamiltonian. The manifold of stationary states is then the span of the Pj. 
The (super) projections V on this manifold and its complement Q, are given 

by 

P(p) = ^P,pP„ Q{p) = J2PjPPk- (11) 

V and Q are orthogonal projections. V satisfies Eq. (7) and Q satisfies 
Eq. (9). 

3 Stationary states 

Let us consider the (super) projections V on the manifold of stationary states 
from a perspective that puts the special classes treated above in a uniform 
context. 

It is a basic property of Lindblad operators that [4] 

Ker£ n Ranged = {0}. (12) 
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Figure 2: The (super) projection V projects on Ker£ along Range £. The 
two spaces are transversal in the Lindblad setting. 

This follows from exp(t£) being a contraction: = implies ex])(tC)p = 
p + tCp and we conclude Cp = 0. 

This allows to define a projection on Ker(£) © Range(£) by- 



Assumption 2 (Gap condition). is an isolated point in the spectrum of 
C and V is given by the Riesz projection 



where the contour encircles but no further points of the spectrum (see 



Remark 1. If is an eigenvalue of finite algebraic multiplicity, then the 
Riesz projection part is for free. The assumption is implied as a whole by the 
strong version of Assumption 1. The assumption guarantees that Ranged is 
a closed subspace. 

The consistency of Eqs. (13) and (14) deserves a discussion. In fact, the 
Riesz projection V always satisfies the first line of Eq. (13) and the validity 
of the second one is the core of the assumption. To see this consider, besides 
of V given by Eq. (14), also Q similarly given in terms of a contour encircling 
the complementary part of the spectrum. Then 



Vp = 



p when p G Ker(£) 
when p G Range(£). 



(13) 




(14) 



Ftg. 3). 



V+Q = l, [P,C] = [Q,'C] = 0, 



(15) 



7 



Figure 3: Since £ is a contraction, its spectrum away from the origin is 
contained in the elhpsoid blob in the left half-plane. The origin is assumed 
to be an isolated point in the spectrum. The circles are the integration 
contour for the Riesz projection V and Q. 



proving the first line. Now assuming the second line of Eq. (13), the eigenspace 
associated to has a trivial Jordan block. This means that the Laurent ex- 
pansion of the resolvent does not have a z~'^ term and hence 

VC= — I dz^- = (16) 

and 

QC = C, C* = C*Q*. (17) 
This places Eqs. (7, 9) into their general natural context. 

Example 1 (Gapless Lindbladians). It may happen that H is gapped and L 
is gapless. For example, let H be a Hamiltonian with a ground state separated 
by a gap from a continuous spectrum. Then the associated Lindbladian Cp = 
—i[H,p\ has the eigenvalue embedded in the continuous spectrum. 

4 Controlled Lindbladians 

The geometric aspects emerge when one turns one's attention to a parame- 
trized family of Lindbladians £<^. We shall call the parameters G 
control and M. the control space. This makes the Hamiltonian i?(0) and the 
coupling to the bath Ta{(j)) functions of the controls. The explicit form of 
these functions is, of course, model specific. 

Assumption 3 (Controlled Lindbladians). 
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(A) The Lindbladian is a hounded (super) operator which is a smooth 
function of the controls 0. 

(B) The gap condition, Assumption 2, holds for all L^. 

4.1 Iso-spectral Lindbladians 

A distinguished family of controlled Lindbladians is the family of iso-spectral 
Lindbladians given by the action of unitaries on H and F: 

i/((/>) = [/(<)6)i/[/*(0), r„(0) = t/(0)r,[/*(0). (18) 

The Lindbladian describing a harmonic oscillator coupled to a thermal bath, 
whose anchoring point is controlled, is an example: 

Example 2 (Controlled oscillator in thermal contact). A Harmonic oscil- 
lator, anchored at the origin and coupled to a heat bath, is described by the 
Lindbladian 

H = a*a, r_ = V7_«, r+ = V7+«*, (7- > 7+ > 0) (19) 

and where y/2a — x + ip. The stationary state of the oscillator is a thermal 
state with (3 = log(7_/7+) [13]. The Harmonic oscillator with controlled 
anchoring point is described by the iso-spectral family with U{4>) — e"'^*^. 
Explicitly 

V2a{<j)) = V2U{(f))a U* (0) = {x - (p) + ip. (20) 

Since the T 's adjust to H the oscillator wants to relax to the thermal state of 
the instantaneous Hamiltonian. 

4.2 Parallel transport 

We shall denote instantaneous stationary states by a. By definition V{a) — a 
(we allow dimT' > 1). By Assumption 3 the projection on the stationary 
states, is a smooth projection on control space and dV is a bounded 
operator valued form. 

The differential of o" = Va gives the identity da = {dV)(j + Vda . Since V 
is a projection V{dV)V — and consequently Q{da) = {dV)a, while V{da) 
remains undetermined. Parallel transport is the requirement, given in two 
equivalent forms, 

V{da) - 0, da^ {dV)a. (21) 
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This evolution of a is indeed interpreted geometrically as parallel transport: 

There is no motion in V. The case dimP = 1 is a special simple case, in 
that there is a unique state cr = (t(0) in the range of It solves Eq. (21) 

without further ado. 

Proposition 1. Under Assumption 3, the form da is trace class. 

Proof. Follows from Eq. (21), the fact that a is trace class and the bound- 
edness of dV. □ 

4.3 Holonomy of parallel transport 

In general, parallel transport, Eq. (21), does not integrate to a function 
on control space Ai unless the curvature vanishes: VdV A dW = (see 
Appendix A). If such a function a = (j(0) exists, it will be called an integral 
of parallel transport. This is, of course, automatic if either = M, or 
dimT' = 1, (see Eq. (5)). 

Parallel transport is consistent with the convex structure of stationary 
states [4]. As a consequence it preserves extremal stationary states. Recall 
that a (stationary) states is called extremal if it can not be written as a convex 
combination of two other (stationary) states. For such extremal states we 
have: 

Proposition 2 (Parallel transport of extremal states). The parallel trans- 
port equation takes extremal stationary state to an eoctremal stationary state. 
If, moreover the manifold of (instantaneous) stationary states is a simplex, 
spanned by a finite number of isolated extremal states (7j{(t>), then the function 
on M 



with pj independent of (j), is an integral of parallel transport. 

Proof. A more general statement has been proved in [4, Proposition 3]. The 
intuition is that states move inside Ker£ as little as possible. In particular 
the boundary should be mapped by parallel transport to the boundary and 



Parallel transport is path independent for two important families of Lind- 
bladians: 

• Lindbladians with a unique stationary state, (Section 2.2). 




(22) 



extremal points to extremal points. 



□ 
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• Dephasing Lindbladians where the isolated extremal states are the one 
dimensional spectral projections (7^(0) = -Pj(0), (Section 2.3). 

Finally we discuss the parallel transport for iso-spectral families. 

Proposition 3. The unitary family o"(0) = U {(j))aU* {(j)) , where Gd(j) = 
iU*dU, is an integral of parallel transport, Eq. (21), if and only if 

V{[G, a]) = 0. (23) 

In particular this is the case when u is an isolated eodremal point. 

Proof. Condition (23) follows by inserting 

da = -i[G, (T]d(t> (24) 

into Eq. (21). The last claim is a consequence of Prop. 2. □ 

Example 3. Condition (23) holds for any iso-spectral family with a unique 
stationary state (and G bounded). In this case V is given by Eq. (5) and 

V([G, a])^aTr([G, a]) = 

by the cyclicity of the trace. 

5 Observables and fluxes 

We denote observables by X. The evolution of observables (in the Heisenberg 
representation) is generated by C*: 

X^dtX-rC{Xl (25) 

where 

c*{x) = i[H,x]+v*{x), v*{x) = ^r;[x,r„] + [r*,x]r„. (26) 

a 

X is itself an observable: We refer to X either as the flux (or rate) of X or 
simply as the flux X. For example, the velocity is the flux of the position 
and the force is the flux of the momentum. 
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Assumption 4. X is not explicitly time dependent, i.e. dfX — 0, and hence 
X^C*{X). 

Fluxes lie in Ranged*. They have the special feature of vanishing expec- 
tation in stationary states. In fact: 

Proposition 4 (No currents). Let a be a trace class stationary state and 
X the flux of the bounded observable X with H and Fq, bounded. Then, 
Ti[Xa) = 0. Conversely, ifTi{Aa) = for any such stationary state, then 
A — C*{X) = X for some bounded observable X. 

Proof. The (super) operator C* acts on the space of bounded operators, this 
being the dual of the space of trace class operators. By Eq. (25) 

Tr(Xa) = Tr(£*(X)a) = Tr(X£(a)) = . (27) 

Conversely when (a bounded) A has vanishing expectation in stationary 
states, then 

= Tt{AV{p)) = Tt{V*{A)p) (28) 
holds for any p. Hence 'P*{A) = and A lies in the range oi £*. □ 

An example of an observable which is not a flux is: 

Example 4 (Hopping on a clock). Consider a quantum particle on a ring 
with p sites, Zp, evolving by the (bounded) Hamiltonian 

H{(j)) = e^'^T + e-^'^T*, {T-^){n) = ^(n- 1). 

p(f) may be interpreted as the magnetic flux threading the ring. The stationary 
states are {n\ aj \m) = p-ie*(""™)2vri/p^ (j = 0, . . . p — 1). The angular velocity 
is the (bounded) operator 

-d^H ^ -i{e''^T - e-''''T*), 

Since —Tr:{d^Haj) — 1€\Tii^'K2 jp — 0) does not vanish in stationary states, 
d^H is not a flux: The angle is not an observable. 

The velocity is the flux of the position operator, which is usually un- 
bounded. It is therefore interesting to examine conditions that would allow 
extending Prop. 4 to unbounded operators. Indeed, the gap condition im- 
plies that the expectation values of fluxes vanish in stationary states even for 
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unbounded X provided Q*{X) is bounded. Indeed, the gap condition allows 
us to use Eq. (17) and replace Eq. (27) by 

Tr(£*(X)(T) = Tr(£*(Q*(X))o-)) = Tr(Q*(X)£(a)) = . (29) 

A more careful discussion of this point is given in Prop. 16 of Appendix B. 
Historically, equilibrium currents in superconductors and mesoscopic systems 
were a matter of debate, see e.g. [11, 12]. 

An example where X is unbounded but Q*{X) is bounded is: 

Example 5 (Taming X). Consider a quantum particle with spin hopping on 
the integer lattice. The Hilbert space is ® and let the Hamiltonian be 

where T is the unit left shift and a and a* are the spin lowering and raising 
operators — (a*)^ = 0, {a, a*} = 1. Since H"^ — \ we can write H as a 
difference of two (infinite dimensional) projections: 

H^P+- P_, 2P± = 1 ± H. 

The position operator is X ®\. It is clearly unbounded. Since [T, X] — T , 
the velocity is the bounded operator 

X = i[H, X®l]^i{T®a-T* ®a*). 

In fact, X'^ = 1. The appropriate version of Eq. (11) says that Q* is given 
by 

Q*{X) = P+XP_ + P_XP+ 

= P+[X,P_]+P_[X,P+] 

^l(-P4X,H]+P4X,H]) 
being a product of bounded operators, it is bounded, even though X is not. 
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5.1 Virtual work for Lindbladians 

The principle of virtual work associates observables with variations of a con- 
trolled Hamiltonian H{(/)) (see Section 4). Our aim here is to formulate a 
corresponding principle for Lindbladians. 

Observe first that £ is a (super) operator, so its variation does not define 
an observable and moreover, the notion of "energy" is ambiguous in Lindblad 
evolutions. The principle of virtual work we formulate is gauge invariant in 
the sense of Section 2.1. 

Theorem 5. The observables given by the joint variation of H and V 

x^sr = SH + zJ2 (r;5r, - 5r;r„) (30) 

a 

are (formally) self-adjoint and free from the ambiguity in H and V, under 4> 
independent gauge transformations. 

Proof For g^ and e independent of (p, the gauge transformation, Eq. (3), 
affects the variation by 

{SH, Sr^) ^(SH-i Y,i9*JT^ - gJTl), 5T^) . 

a 

This leaves X^ invariant. The same applies to the transformation (4). □ 

Remark 2 (A second family). A second family of observables that are gauge 
invariant is 5 {Y.a f^a, T*]). 

The observables X^ extend the notion of the principle of virtual work to 
the Lindbladian setting. The physical interpretation of X^ is often suggested 
by dimensional analysis. However, their interpretation depends on the precise 
choice of controls and is model dependent. 

Example 6 (Controlled oscillator in thermal contact: Example 2 continued). 

Shifting the anchoring point of the oscillator gives 

SH - -"^^H = -xs<p, t{r*jr± - 6r*^r±) = ^^^-^54, = ±^^p5(t>. 

—X is the spring force, while — 7_p gives the friction force due to the cold con- 
tact and 7-i-p the gain from the hot contact. The total force is the momentum 
flux 

p = - (7_ - 7+)p, (7_ > 7+). (31) 
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In the example, the principle of virtual work gives a flux. This is not a 
coincidence. For iso-spectral Lindbladians, Eq. (18), virtual work is a flux. 
More precisely, let denote the (local) infinitesimal generators 

G5(t) = G^5(t)i' = iU*5U (32) 

(summation implied). The variations are: 

5H = i[H, 5r, Sr^ = i[r„ 5r- (33) 

Theorem 6 (Virtual work and fluxes). For iso-spectral families of Lindbla- 
dians generated by G^, the observables associated with the principle of virtual 
work, Eq. (30), are the fluxes of the generators G^ : 

c*{G^)5r ^SH + tJ2 (K^^e. - 5r;r,) . (34) 

a 

In particular we have Noether's theorem in the form: IfSU is a symmetry, in 
the sense that the r.h.s. vanishes, then its generator is a conserved quantity. 

Proof. Transparent. □ 



5.2 Currents 

Just as there are several notions of force in a damped oscillator, there are 
several notions of currents in an open system. In an open system charge 
need not be conserved as charge may migrate between the system and the 
bath. One can distinguish three notions of currents, related by one relation: 
The rate of charge in a subsystem; the rate of charge transfer between a 
subsystems and its complement, and the rate of charge transferred from the 
subsystem to the bath. 

More precisely, consider a subsystem associated with a domain Q. Let 
Qo, denote the observable associated with the charge in VL. The first notion 
of a current is the flux of Qn, namely, 

= C*{Q^) = i[H, Qn] + V*{Qn) (35) 

Qci, being a flux, is gauge invariant in the sense that it depends only on C 
and not on its partition into H and T. It is, however, non-local in general 
since Fq,, in contrast with H, need not be local. 
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We assume that H is charge conserving, i.e. [if, Qnun^] — 0- This allows 
to define the current from the subsystem Q to its complement by 

Iau = t[H,Qn]. (36) 

Since H is local, is localized near the boundary, dVt^ between the subsystem 
and its complement. However, it depends on the partitioning of C into H 
and Fq. 

The remaining term gives the rate of dissipated charge, Sq. By definition, 
the three currents are related by 

Qa = Ida + Sa. (37) 

The three currents have different characters and are measured by different in- 
struments. Qq, is measured by an electrometer while Iqq, can be measured by 
an ammeter that monitors the flow at the boundary between the subsystems. 
The last term has been called dissipative current in [8, 16]. 

Example 7. Consider Fermions hopping on a one dimensional lattice which 
can also tunnel in and out of a bath. The Lindbladian has 

j 

with ttj the usual Fermion annihilation operators for site j. The charge in 
the left semi-infinite box is 

j<0 

and the currents in Eq. (37) are 

Iql = i{a\ao - aoOi), -S'l = 2 ^(7+0^0* - 7_a*aj) : 

Iql is localized at the boundary of the box, whereas Sl is not. 

Dissipating currents can arise also when the Lindbladian is charge con- 
serving as we discuss below. 
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5.3 Currents in a magnetic field 

The action of a magnetic field on cfiarged particles endows the dynamics 
with chirality. This has interesting consequences for currents. Consider the 
Lindbladian describing a charged particle in the plane under the influence of 
a constant magnetic field coupled to a heat bath. The Hamiltonian is the 
Landau Hamiltonian^: 

H^D*D, D ^ -idi + 82 + Bx2 = vi + iv2, (vi^^vl), (38) 

and the (Markovian) thermal bath, is described by (cf. Example 2)) 

r_ = ^D, r+ = ^D*, (7_>7+>0). (39) 

We shall call the generator of the corresponding evolution a thermal Landau 
Lindbladian. 

Proposition 7. The (total) current density of the Landau Lindbladian of 
Eqs. (38, 39) is 

Jm(^o) = {p{xo),Vfj} - (7+ + 7_)a^p(xo) + 2(7_ - 7+)p(a;o)£^^'t;^. (40) 

The charge density is p{xq) — 5{- — Xq) and is the completely anti- 
symmetric (Levi-Civita) tensor. The (total) current satisfies charge conser- 
vation: 

dtp = -d,j,. (41) 

Before proving the statement, let us comment about its content. The 
Hamiltonian current is proportional and parallel to the velocity 2f^. The 
dissipative current has a (non-chiral) diffusive term proportional to the gra- 
dient of the density and a further chiral term. The dissipative currents can 
be interpreted in terms of Brownian motion (see below). 

Proof. The dissipative terms of the Lindbladian are 

vi{x) = r^[x,r±] + [r*^,x]r±. 

For a function X — f{xi, X2) of position we have 

i[H,f]^{v„dJ}, Pi(/)=7±A/T27±(9^/)£^,^.. (42) 

^Strictly speaking, this example lies outside the scope of our framework since it in- 
volves unbounded operators. However, it is simple enough that one can check that formal 
manipulations are indeed justified. 
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Eq. (35) then gives the dual form of the statements of the proposition, namely, 
dtf = J jix{xo)d^f{xo)d'^Xo, f = J p{xo)f{xQ)d^XQ. 



□ 



Stochastic interpretation 



The dissipative currents admit an interpretation in terms of a (classical) 
stochastic process. To see this note first that for functions X — f{vn) of 
either velocity {/i — 1,2) 

T^*±(f) = l±B'f" ± 2^^Bfv, (43) 

which can be read as if originating from a (exciting or damping) Langevin 
equation 

dVfj, = ±2-f±Bv^dt + Bdh^^t, 
where is a Brownian motion with zero drift and variance 

^{dh^^tdhy^t) = 2j±Sf,i,dt. 

In fact, expanding E(/(t'^ + dv^)) to first order in dt and to second order in 
dbt yields that expression. 

To derive the Langevin equation for dx we first note that the guiding 
center (ri,r2), 

r,-x,+ '-^v. (44) 

satisfies [r^, v^,] — and thus is a constant of motion for the Lindbladian, 
— 0. Insisting on being a constant of motion, we have 

dx^ = -^dvy = e^u{T'2l±Vudt - dK^t)- (45) 

JD 

In view of ¥.{e^ydh^^t£^iyidhyi ^t) = '2l±^iJ.fj,'dt this is the Langevin equation 
corresponding to Eq. (42). (Beware: dx^ ^ v^dt.) 

Wc can now combine C*{r^) = with Theorem 6 to conclude 

Proposition 8 (Iso-spectral Landau Lindbladians). The velocity, x, is the 
(negative) virtual work associated with the iso-spectral family of Landau Lind- 
bladians {H{(j)),T±{(l))} generated by 

= B-'e^,v,. (46) 
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The unitary acts on wave functions by 

(C/^V)(^i, ^2) = H^i + h/B, X2 - <l>rlB)e"l>^^\ (47) 

The physical interpretation emerges by noting that 

U^v^U; =v^- 0^. (48) 

Since (f)^ appears in H hke a pure gauge field, its variation in time, —0 is a 
constant electric field that drives the system. Alternatively, the proposition 
may be viewed as a manifestation of gauge and translation covariance, in the 
sense that —id^ and appear in the Lindbladian only through the minimal 
coupling expression v^. The variations in Eq. (34) generated by — G^n and 
are then the same by 

in fact both sides equal S^i,. 



6 Adiabatic Response 

We are interested in adiabatically changing controls; (j) — (p{s) where s — et 
is the slow time. The evolution equation for the state p is 

^1 = C,p. (49) 

with initial state that is an instantaneous equilibrium state. A key feature 
of adiabatic theory is that the evolution of p is slaved to the evolution of a. 
We borrow from [4]: 

Proposition 9 (Adiabatic evolution). Under Assumption 3 the solution of 
Eq. (49) with initial condition the stationary state a{0) is 

M \ i \ , f 0(6'^) z/dimT'^l 
(T'P)(^)--(^) + [oU z/dim7^>2; ^^^^ 

and 

{Qp){s) ^ eC-'a{s) + 0{e% (51) 

where a{s) is the corresponding integral of parallel transport. 

C~^{a) is well defined and bounded since a e Range £. This follows 
from parallel transport a — Q& and the definition of Q as the projection on 
Range C 
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6.1 The response of unique stationary states 



We are interested in the response of the observable X of an adiabatically 
driven system. The case of a unique stationary state is simpler than the 
general case and we treat it first. 

Proposition 10 (Response coefficients). Suppose Assumption 3 holds and 
the stationary state is unique. Let X he a hounded ohservable and p a solution 
of the adiahatic Lindhlad evolution, Eq. (49), with initial state a normalized 
stationary state cr(0). Then, the response at slow time s is memory-less and 
is given hy 

Tr{Xp{s)) = Tr{Xa{<f>)) + sF^ + ^(s'), 

(summation implied) with cr{(f)) the instantaneous stationary state and (f) = 
(t){s) e M.. The response coefficients 

F.(0) = TY{XC-^\d,a{(t>))) (52) 

are functions on control space M.. 

Proof. This is a direct consequence of the adiabatic expansion, Eq. (50). d^a 
is trace class by Prop. 1. □ 

The first term Ti{Xa) is of 0(1), and describes the persistent response., 
a property of the stationary state. The second term is the driven response 
which is proportional to the driving £0, the (unsealed) velocity of the controls. 

One is often interested in situations where F[<^) is constant on M.. This 
feature depends on additional structure (e.g. thermodynamic limit, disorder 
[9, 6, 1]). Observe that the expression for F(0) involves inverting £, an 
operator with a non-trivial kernel, and so is not completely elementary. 

Remark 3. The formula for the response coefficient F^, can he cast in a way 
that is formally reminiscent of Kubo 's formula: 

POO 

F^{<P) = - lim / dm((e(^^"^)*X) (a^cr(0))). 
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7 Response of fluxes 



When the manifold of stationary states is multidimensional, the persistent 
response has memory and F can not be viewed anymore as functions on 
control space M. (see Section 4.3). However, in the case of observable which 
are fluxes several simplifications occur: There is no persistent response and 
the formula for F„ simplifies and becomes elementary. If, in addition, the 
extremal stationary states are isolated (Section 4.3) then, in addition, F 
defines a function on M.. 

By definition, a fiux (which is not explicitly time dependent. Assump- 
tion 4) can be written as 

sX = C*{X) = C*{Q*X), (x^^y (53) 

where the replacement X by Q*{X) relies on Eq. (17) and is only of interest in 
the infinite dimensional case where X is unbounded while Q*{X) is bounded. 

Theorem 11 (Response of fluxes). Suppose Assumptions 3 hold; that (t(0) is 
an integral of parallel transport; the flux X and Q*{X) are bounded operators. 
Then, to leading order, the response is memory-less, linear in the driving and 
given by 

Tr{Xp)(s) = FM + 0(s), = ^) 

(summation implied) and where the response coefficient 

FM = Tr(Q*(X)a,(7) 

is a function on Ai. 
Proof. From Eqs. (15, 51) we have 
Tr{Xp) = e-^Tr{/:*{Q*X)p) = Tr (£* (Q*X)£-V) ^ Tr {Q*{X)d,a) 0^ 

□ 

The observation that one can sometimes avoid computing Green func- 
tions in linear response is at the heart of the TKNN formula for the Hall 
conductance [24]. 
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7.1 Geometric magnetism for iso-spectral families 

The response coefficients of an iso-spectral family generated by are nat- 
urally organized as a matrix F^j,, relating the response of the flux of to 
the driving (ff . The analog of the formula in Prop. 10 is 

IV(G^p)(s) = F^,(0) 0^ + 0(5). 

Combining Theorem 11 and Eq. (24) we get for the response matrix 

i^M. = IV(Q*(G'^)9.a) = -iTx{Q*{G^)[G,,a]). (54) 

We are now ready to state our main result: 

Theorem 12 (Geometric response). Suppose G^ are hounded and a{(f)) — 
U {(f))(TU* {(f)) is an integral of parallel transport. Then the response matrix is 
antisymmetric and given by 

F^,^-iTr([G^,G,]a). (55) 

//, moreover, a is a projection P then F is the adiabatic curvature of the 
bundle (see Appendix A): 

F^,^iTr{P4d^P,d,P]). (56) 

For unitary evolutions, the first part of the theorem reduces to a (special 
case of) result of Berry and Robbins [10], who coined the term geometric 
magnetism for the anti-symmetric part of F. 

The second part of Theorem 12 extends the geometric interpretation of 
response matrix from the unitary case [7] to open systems. The conditions 
in the theorem are satisfied for Lindbladians representing relaxation to the 
ground state and dephasing Lindbladians whose initial state is a spectral 
projection. 

Proof. The conditions have been set so that the formal manipulations are 
justified 

F^, = Tr{G^d,a) = -iTr{G^[G,,a]) = -tTT{[G^,G,]a), 
where in the second equahty we used Eq. (24). 
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For the second part observe that the equation —i[G, P] — P imphes 
P±GP = iP±P and PGP± = -iPP±. 

Hence 

F^, = -tTT{[G,.^ G,]P) - -iTr{PG,P^G,P - PG.P^G^P) 
= -iTi{d^PP^d,P - d^PP^d^P) 
= zTr(P^[a^P,a,P]), 

where the first hne is a readily checked identity. □ 

Remark 4. More details about the dephasing case are in Appendix A. In 
particular, we give a formula for response when (7{(t>) is not the integral of 
parallel transport. 

In the case that [G^, Gi] is proportional to the identity, the transport co- 
efficients are purely geometric and independent of the dynamics. An example 
of this kind is the Hall conductance on a torus^. 

Example 8 (Hall conductance). Consider the iso- spectral family of Landau 
Lmdhladians, {H ,V ±{<p)} of Prop. 8, with B = 2n so [D,D*] = Att. By 
considering the Hilbert space of quasi-periodic functions 

i>{xi, X2) = i>{xi + 1, X2) = e^""^^ ^(xi, X2 + 1) 

one gets an iso-spectral family of Landau Lindbladians on the unit torus. 
(The space is invariant under the action of U^). Their stationary state is 
the therm,al state with An 13 = log7_/7+. By Eq. (48) changing the con- 
trols {(f)i,(f)2) iTiay be interpreted as changing the two Aharonov-B ohm fluxes 
threading the two "holes" of the torus, Fig. 4- The response coefficients have 
the meaning of conductances. By the first part of Theorem 12 

F^, = 27r£^, (57) 

independent of'y±. ( The insensitivity of the Hall conductance to temperature 
is a pathology of the model.) In the units we use 27r is the quantum unit of 
conductance. 

An illuminating example of Theorem 12 in the case when [G^, Gi,] is not 
proportional to the identity was recently given by Read and Rezayi [20]. 
When the G^j arc the generators of shears, their commutator is a rotation, 
hence relating the Hall viscosity to the angular momentum. 
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Figure 4: A two dimensional torus threaded by two Aharonov Bohm flux 
tubes carrying fluxes 0i^2 that serve as controls. 

7.2 Friction and dissipation 

The fact that F^i, of Eq. (55) is anti-symmetric does not imply the absence 
of dissipation. It only says that looking at the response of fluxes is not ap- 
propriate for the study of dissipation. To explain this statement consider the 
dissipation associated with the dragging of the anchoring point of a (damped) 
oscillator coupled to a heat bath at velocity 0. The response coefficient relat- 
ing force to velocity is friction. As there are three forces in the problem — the 
momentum rate, the force on the anchoring point, and the friction force — 
there are also three friction coefficients. The friction coefficient associated 
with the momentum rate vanishes, but the others do not. 

Example 9 (Friction: Example 6 continued). The (unbounded) generator 
of shifts is, p = i{a* — a)/\/2. With a a thermal state of the oscillator, 
d^a = —i[p,a] is trace class and Eq. (55) applies with = Gy = p. The 
friction coefficient vanishes, as it must by anti- symmetry. 

The momentum rate vanishes because it can not disentangle the heat lost 
to the bath from the mechanical work done by the anchoring point. To study 
dissipation it is not enough to look at the response coefficients of fluxes, nor 
is it enough to examine the energy of the system. 

Indeed, the energy of the (small) system, in the adiabatic limit, is 

E = Tt{Hp) = TT{Ha) + 0{e) (58) 
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by Eq. (50). In particular, for an iso-spectral family the energy is constant 
(to leading order) when H and a undergo the same unitary transformation, 
as is the case in the example of the damped oscillator. The energy does not 
reveal the dissipation. 

To reveal the dissipation one needs to look at the breakup of the energy 
to work and heat. The variation of the energy 

5 Tr{Hp) = Tr{H 6p) + Tr{6H p) (59) 

expresses the first law of thermodynamics [22] 

5E = 6W + 5Q = {Ti{ad^H) + Ti{H d^a))6(t) + 0{e). 

To compute the friction one needs to study the expectation of the spring force 
—X rather than the momentum flux C*{p). (More generally, d^j,H rather than 
the flux£*(G'^) of Eq. (34).) 

In general, the computation of Ti^pd^H) is complicated for two reasons: 
First, one needs to evaluate C~^. Second, in the case that the ground state 
is non-unique, it also needs the explicit expression for the 0{£) term in the 
adiabatic expansion, Eq. (50), which are history dependent. For dephasing 
Lindbladian such a computation is given in [5]. We shall not pursue this 
direction here. 

As a sanity check, let us derive the first law of thermodynamics using the 
tools of the previous sections. Since H is explicitly time-dependent Assump- 
tion 4 does not hold for its flux is now made of two terms: 

H^{d^H)r + r,*{H). (60) 

Substituting in Theorems 10, 11 indeed reproduces the flrst law: 

^^{TT{ad,H) + Tr{Hd,a))[^)+0{e^). (61) 

8 Concluding remarks 

We have derived a simple and general formulas for the adiabatic response 
coefficients for observable of the form X — jC*(X). In the case of iso-spectral 
families of Lindbladians, the response matrix is determined by geometry 
and is purely anti-symmetric. We find a range of circumstances where the 
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response coefficients are given by tfie adiabatic curvature of tfie associated 
stationary projections. It will be interesting to extend the theory to models 
of extended systems with (non-interacting) fermions. 

Acknowledgments. JEA is supported by the ISF and the fund for 
promotion of research at the Technion. MF was supported by UNESCO and 
ISF. We thank M. Porta for useful discussions. 

A Geometry of projections 

Consider continuous orthogonal projections Pj{4>) with Yl-^ji^) ~ The 
superprojection that takes p to Range P = Span{Pj} is, Eq. (11), 

We are going to describe parallel transport inside Range P [18]. 

For a given path Pj{t), parallel transport ii — PjPjU maps vectors u{0) in 
the range of -P,(0) to vectors u{t) in that of Pj{t). That map U{t) is unitary 
and generated by 

K:^iU*U^Y.'^^^ ^J='PjPr (62) 
j 

In fact K* = K, since A* = -iPjPj = -iPj{l - Pj) = -iPj + Aj, and, for 
U so defined, u{t) = U{t)u{0) satisfies 

ii — —iKu — PjPjU, 
j 

as required. And for a{t) — U{t)aU*{t) the parallel transport equation (21), 

V{t)a{t) = 0, holds true. 

When dimP^- = 1, the parallel transport is manifestly path independent. 
In general, tliis is determined by the standard condition of vanishing curva- 
ture: 

Proposition 13. Let A = {dV)V be an operator valued 1-form. The differ- 
ential equation 

da — Aa 
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admits a (locally path independent) solution a if and only if the curvature 
vanishes 

n = 0, n = -iPdV AdVP. (63) 

It implies the following criterion for the case of parallel transport of pro- 
jections. 

Proposition 14 (Adiabatic curvature). The parallel transport constructed 
above is locally path independent if and only if the adiabatic curvature 

j 

commutes with all elements in Range 7^. 

Proof. Parallel transport of a vector Uj in the range of Pj along an infinites- 
imal square d(l)^d(f)^ maps 

% ~^ Uj — iR^jii,Uj<\(t)'^d(f)^ -\- o{d(f)^). 

The associated adjoint transformation maps the state a = PjcPj as 

a a - i[R^iy, ajdcji^dcj)" + o{d(i)^). 

This allows to read off the curvature TZfj,,, of P seen in Eq. (63): By TZ — TZP 
we have 

PfiuP = [R^,u,Pp]- 

Hence the criterion of vanishing curvature states that R commutes with all 
elements in Range 7^. 

Computation gives the commutator 

[K„K] = Y^P,^[d,P„d^P,] (64) 
j 

as the sum of adiabatic curvatures of all the spectral projections. Since 

j 

one finds 

j 

□ 
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For an iso-spectral family of projections 

Pj{(j)) = exp{-iG(j))Pj{0) exp{iG<j)) 
the generator of parallel transport, Eq. (62), is 

K^G-J2 P3GP3 ^G-V*{G), 
j 

since V* = V hy Eq. (11). While it does not coincide with G it differs from 
it only inside Range V 

Q*{G) = Q*iK). (65) 

When rankPy > 1 and dimA^ > 1 the parallel transport can not be 
integrated in general. And the response coefficients are not functions on the 
manifold. 

Theorem 15. Suppose G^ are hounded and C is a dephasing Lindhladian. 
Then the response associated to the driving path (f){s) and flux G^ depends 
only on the integral of parallel transport a{(f)) and the derivative 5(f) at the 
end point, 

Tr{C*{G^)p{s))^F^,d,cf>{s), 

where 
Proof. 

F^, = Tr{Q*{G^)d,a) 

= i Tr{Q*{K^)[K,, a]) = iTr{[K^, K,]a), 

where the second line express parallel transport and uses Eq. (65). The last 
equality is by V{[K^,a]) = 0, which characterizes parallel transport and by 
the way restates Eq. (23). □ 

Example 10 (Taming X: Example 5 continued). Consider a family of 
Hamiltonians generated by a momentum shift 

^-i^x jj^i<t,x ^ ^ ^ ^ ® a* = P+(0) - P_(0), 

where X is the position operator. The generator of the parallel transport is 
K = ^(P+P+ + P_P_) ^-'-[H,H]^ a*a - ^. 
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Although K ^ X , their difference commutes with the Hamiltonian. Further- 
more K intertwines P±, 



which is equivalent to the statement that K generates no motion inside Range P±, 
P±KP± = 0. 

B Currents and unbounded observables 

We discuss the precise meaning of Eq. (29) when X is unbounded. We still 
assume that H and r„ are bounded, while X — X* need not be. Yet, the 

commutators [i?, X] and [r„,X], defined as quadratic forms on the domain 
D{X) of X, arc assumed bounded, and V^D{X) C D{X). Then X = C*{X) 
is a bounded operator by natural interpretation of Eq. (25) in the sense of 
quadratic forms. 

Proposition 16. Under the stated conditions, Ti{Xa) = for any (trace 
class) stationary state a. Moreover, Q*{X) is well-defined as a bounded 
operator. It is given as a strong limit, Q*{X) — s— lim„_^oo Q*{^n)) by means 
of any sequence of bounded approximants Xn with Xnf — >■ X(p, ((f e D{X) ); 
finally C*{X) ^ a {Q*{X)). 

Proof. There exist sequences X^ as stated, e.g. X^ — X/{1 + n~^X'^). The 
assumption states that the bounded operator [H, X] is characterized by the 
property 

{cp, [H,X],I;) = {Hip,X^) - {Xip,H^), {cp,^ e D{X)) . (66) 

Hence [H,Xn] [H,X] (weakly), and similarly for [ra,X]. Thus C*{Xn) 
C*{X) (weakly). Using that An ^ A (weakly) and B trace class imply 
Tr{AnB) Tr{AB), we conclude 



by Eq. (27). Moreover, by Eq. (17) we have C*Q*{Xn) C*{X). We notice 
that C* is weakly continuous, and so are {C* — z)~^ and the inverse of C* on 
Range £* — Range Q*, i.e. 



KP^ 



P_K. 



Tr{C*{X)a) = lim Tr{C*{Xn)a) = 



(67) 




(68) 
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in the notation of Eq. (14). As a result Q*{Xn) is weakly convergent to a 
limit denoted Q*{X), and the result follows. □ 
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